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FLUID MECHANICS OF AXlSYMMETRlC WAKES 
BEHIND BODIES IN HYPERSONIC FLOW 

By Lu Ting *nd Paul A. Libby 

SUMMARY 

In this rcportTfti analysis of the fluid mechanics of wakes behind 
bodies in hypersonic flow is presented* The models used for idealising both 
the flow field and the chemistry of air are discussed first. The fluid 
mechanics of the wake itself is then shown to be describable to terms of 
transformed coordinates which are applicable to both laminar and turbulent 
wakes. The transformation to the physical plane requires description of 
the transport properties of the gas; these are therefore discussed next. 

In this connection a rational means for estimating the turbulent transport 
properties is described. Finally, there are presented the results of 
numerical analysis of four cases which appear to lead to the highest 
concentration of electrons in the wake.and of the initial conditions for a 


boundary layer induced wake 
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SY MBOLS 

a amplitude o* the velocity dittribution as appears m Equation (14) 

b amplitude of mass fraction distribution 

d distance oi the cross section of the wake from the body 

0 base diameter of the body 

f(h) transformed stream function 

h statu enthalpy; h » 0 at J K 

h* dimensi'ideiss quantities with respect to free stream conditions 

L total entnalpy including mean kinetic energy 

Le Lewis number 

M Mach number 

n electron density 

p pressure 

Pr Prandtl number 

Re Fevnolds number 

R 0 nose radius 

s transformed coordinate along body surface after Lee's theory 

T temperature 

u,v velocity component* respectively in the axial and in 

the radial directions 


free stream velocity 
molecular weight 
cylindrical coordinates 

dimensionless quantities with respect to the radius of the 
initial section of tfje incompressible wake ' 






GCNERAt APMJ1D SCIENCE LABORATORIES, INC 


V 
x 

V 

6 

6* 

« 

n 

Q 


A 

I* 

“t 

t 

P 

iL_ 


**•* j 

mass fraction 

incompressible radial coordinate as defined in Equation (8a) 

ratio of specific heats, -F 

1 v 

radium of the compressible wake 

radius of the incompressible wake as defined in Equation (8e) 
absolute incompressible eddy viscosity as defined in Equation W) 
dimensionless variable as defined m Equation (8b) 
momentum thickness 
l - a 

absolute viscosity 

absolute eddy viscosity as defined m Equation U8) 
dimensionless variable as defined m Equation (J2) 
density 

—stT earn Junction---— 


Subscripts 

40 free st"eam conditions 

e quantities at the outer edge of the boundary layer 

f frozen flow 

1 species i 

j atomss of species j 

it molecules of species k 

° conditions on the plane x » 0 

» stagnation conditions 
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Su per gc ri p t»» 

( ) f time averA^cs of fluctuating terms 

* transformed quantities rtfeVred to turbulent incompressible wake 
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1 INTRODUCTION 

The flow field behind a body in hypersonic flow can be divided into a 
region close iu the body and a far field remote from the body. In the 
former region, the pressure is significantly different from its ambient 
value and is non •uniform in both the radial and axial directions. Moreover, 
in this region the viscous effects associated with base mixing play an 
important role in determining the flow field. The extent of this region 
close to the body can be estimated by application of the blast theory of 
S. C, Lin (Reference l)i if it is required that p-l, l p UJ , with 
v - 1.4, then k/D^ 0.0665 M,o { p/p w , * t*. -40i] l ^46; if p/p^-l then 
x/D54. For the purpose of the wake these values of x^D correspond to 
’•close' 1 to the body# 

The present analysis is concerned with the far field where viscous 
stresses predominate. The static pressure m this region is taken to be 
constant and equal to its ambient value. Therefore, the wake is being 
considered here as a generalization oi the classical wake behind bodies* in 
low speed flow. * In contrast to the low speed case, however, the hyper¬ 
sonic wake can involve changes across the viscous region of stagnation 
enthalpy, composition, and density i - well as of velocity* 

There appear to be two types oi wake in the case of hypersonic flow. 
For blunt bodies the alterations in flow and gas properties from their free 


*TM» problem is reviewed in Reference L 
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stream values are associated almost entirely with the bow shock wave 
and to a lesser degree with the shock wave system near the base of the 
body. The miss flow influenced by the^e shock waves is much greater than 
that involved in the boundary layer and base mixing* Thu* the lateral 
extent of the portion of the wake caused by the boundary layer and base 
mixing is relatively small. The wake in this case can be idealised as a 
constant energ> region with initial flow properties determined by the shock 
wave system. For slender bodies, on the contrary, the wake can be cop* 
side red to be entirely due to the viscous effects associated with the boundary 
layer and base mixing. Then the energy is non-uniform across the walei* 

1 he classification of hypersonic wakes then involves both shock-induced 
and b iundarv layer-induced wakes. In addition, it is possible to consid r 
two different forms of transport or responding to laminar and turbulent flow< 
Finally, assumptions tmiicerning The chemi cal and thermodynamic belMtvJop 
of the gits introduce an additional parameter, A realistic and complete 
description of the behavior of air is complex* In most analyses of fluid 
mechanical problems, it is customary to consider initially certain limiting 
cases of gas behavior. In the analysts described here, the composition of 
air is idealized into molecules, atoms and ions with the ions treated %s 
having a molecular weight equal to that of the molecules. The electrons ^ 

are taken to be equal in number to the number of ion, and to be negligible;T„^ 

v - v 

from the fluid mechanical point of view. Finally, the ion'COae#eEf|Alesi#^*^ 
while of interest, is assumed small so that the flow properties in tke dyl|b 
are determined by the molecules and atoms* ' ‘;^*5§BS 
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With this tertiary mixture as a model, there are various limiting 
cases of chemical and thermodynamical behavior to be considered. The 
usual cases of complete equilibrium and completely frozen composition 
are supplemented by several others* For example, the concentration of 
molecule* and atoms can Le considered to be frozen; from the flow analysis, 
a "state* 1 defined bv a pressure and temperature can be found. The con¬ 
centration of ions in this state can be assumed to be given by equilibrium 
air in that state Such a case corresponds to frozen dissociation and 
equilibrium ionization, A fourth limit ciearty corresponds to equilibrium 
dissociation and frozen icaization. 

Thus the two forms of wake, the two transport mechanisms and the 
four limiting cases of gas behavior result in 16 different types wf wake 
which can be analyzed in terms of the idealizations discussed here The 
analysis of the wake which is presented below is applicable to all of these 
cases; it is baaed on the application of integral methods to wake flow due 
to Bloom and Steiger (Reference <). There result relative!', simple 
expressions for the i‘U«w properties in the w'ake in term* of transformed 
coordinate*, of the gas properties along the axis of the wake, and of 
conditions at the origin of the wake. The application of these expressions 
to a particular wake flow', therefore, involves the determination of the gas 
properties and of the initial conditions for the wake and then the evaluation 
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uf the inverse tran*>formation to obtain the ilo^ parameter* m the physical 
plane 

Primary emphasis has been devoted here b* those ca.*e* which would 
appear to lead to the highest * am enf ration of electrons in the wake . Ptr~ 
lim:nar\ analy-.#-* ;ntH uteci that the sho .jv * induced ^<*kr with complete 
equiiibirum and a; th completely frozen compos.t <>n were thu.» to be 
emphasized, therefore, attention k-. directed to four tases which so 
arise* The theoreim.il of th* bound nr. ln>« r *r.n teed wake is also 

presented «&ioi:*; a w:th a mimerirai mudv-.s ot tin .nitt r.l condition* in ».-ch 

a wake* in th** ;».»». **w 1 section, the initial \ona;;iuns for the shacknucuced 
wake are d; ,u^s<*d 

i he author- -re pirn ed to acknowledge the helpful suggestion# of 
Dr. A11 tom o Kerr; at a; i>r* M TD Bioum * oncer muc the ana-KsD* presented 
here »w;d the v of Mr* Pa us Bart»nt; ,n the numerical ana.lv 
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11 INIT1A1 CONDI HON S FOR THE S HOC K*IN DUCED WAK E 

The major contribut «*?; to the wake behind a blunt body in hypersonic 
flow arises from the bow* shock* It is therefore amumed that the bow shock 
shape is given and that the initial conditions for the wake are derivsi In 
accordance with the idealizations described above, these initial conditions 
involve the radial diaiributionu of the axial velocity, Of the i concentration# 
of molecules, atoms and ions, and of a state variable, for example, tem¬ 
perature, I he statu pressure throughout the wake is assumed constant 
and known. 

From the shock shape and the tables of Feldman (Reference ?}, the 
state of the mixture at each point behind the shock can be found for a given 
flight velocity and altitude. These mixture properties can be usea in con¬ 
junction with the equilibrium compositions tabulated by Logan (Reference 5) 
to determine completely the state of the air behind the shock. An assump¬ 
tion must be made concerning the thermodynamic and chemical process 
experienced by the gas in going from its state behind the shock to the 
ambient pressure p^. For the calculations herein reported, two isentropic 
processes were considered cor responding to complete equilibrium and 
completely frozen flow* With these assumptions, the state of the flow 
(velocity, composition and temperature) where the static pressure is p^ 
can be determined for each streamline originating at the shockwave. There 
remains only the determination of the radial positions of the streamlines. 
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Thi?s i am be clone by a sir earn tub* method uf anaivfti* involving application 
ol the cunt lit) on oi mas* conservation starting from the axis of symmetry* 
For the v of frozen composition the tables oi fhalf el al ;n Reference 6 
were used in * onjunction witu the comp*billon <Uta from Reference % 

I he bodv considered lor name r;cal or.alysjs according to this method 
i» the If) configuration with «* base ratline of one foot assumed. There was 
available a sohlieien photograph of the sbuck hape at a Mach number of b* Q 
(Reference ?). B\ application of the Mach number independence principle 
this same shape was osed to determine the initial value* for the wake 
corresponding to a velocity of 2 j,2O0 'ft/set at an altitude of £00* 000 feet* 
i he results for the initial distributions according to the two aerodynamic 
processes are mown m Figures I to 4, ]t should be noted that it the radial 
distribution# lor a base radius oi the body rqua! to 2* n feet are desired, 
then o u end the radial scale associated with a and r must be multiplied by 
2* 5: that is, the w ake is 2. * umes greater in diameter. 



*««i II 
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111 ANALYS IS OF FLOW IN THE SHOCK- INPUCFD WAKE 

We assume that the wake begins at a plane x * 0 which is at distance 
d 0 behind the body. On the initial plane x 0, all the flow properties are 
prescribed, for example, as shown in Figures l to 4* The problem is to 
study how the flow in the wake approaches the uniform flow outside the wake 
for the esse of shock-induced wake# 

Since the length of wake much larger than the diameter of the cross- 
section, the usual boundary layer assumption is applicable. The pressure, 

p, across the wake is constant and is equal to that of the surrounding uniform 

■■ i 

flow, i. e* , 


P * Poo ( 1 > 

With the assumption that the Prandti number Pr and the Lewis number 

Le are equal to unity, the energy equation for shock-induced wake is fulfilled 


by the Crocco integral 


h ♦ 


2 h« - const, 
s e 


U\ 


The assumptions concerning the Prandti and Lewis numbers are not necessary 
when we employ the integral method in which the differential equations are 
replaced by the integral equations or the conservation laws across the entire 
wake layer. This argument will be substantiated in the treatment of froaen 


The continuity and momentum equations for the wake flow are 


flow 
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and 


pur 


dg 

dx 


+ 


h\ 

PVT - 


d 

dr 


du 


(pr 


dr 


(4) 


These equation* are applicable for laminar flow and also for turbulent flow 
if eddy viscosity in introduced ui place of the visc osity coefficient P, In the 
following analysis it should be understood that the symbol p denotes either 
the laminar or turbulent transport coefficient. When the integral method is 
employed, these two differential equation* for u and v are replaced by the 
integral relationship for u: 


d_ 

dx 


I 


Murlu^- u) dr * 0 


o 


where 6(x) : the radius of the wake regmu. This equation states that the 
total momentum is conserved in the wake t eg:on and that the momentum 


thickness"R“is" a constant '(RefErencel ITahcT <T. Thus 

6_ 


g - - 


1 


1 


P(» U CD 6 0 


pur(u OQ -u)dr- --- ^ ^ 

h ’lX> U cu o t 


1 . 

> \ if Ur < U OD 


-3 


dr* const 9 


tx*o 


{$> 

where 6(x) is the radius, of the wake and 5^ is the associated incompressible 
wake radius at x - 0 as defined in Equation (Sc). 

The boundary conditions for u are 

$xi 
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Ju o u d* u 

and at r - 0 , ^ - 0 and pu - d pr 


ar 


dr £ 


( 7 a . b ) 


i he last boundary condition is obtained from Equation (4) {tf. Reference 8), 

The compressible wake is transiormed into an associated u\e ompressibl^ 
wake by the transformation (Reference 8) 


- \ I. u r* 

J 


which in turn is transformed to the dimensionless variable*?, 

n • * 6 f 


(8a) 


b* 


where 


P .(&')**{ pur* 

<U 'i 


(8b) 

(8c) 


lUi.tions (f>), (»>) f *5 then be<.oit;< 


0. . ( . 

° • 6 • ' l .fl-ildi) 


6* 


J 


(V) 


- i 


a - l. ? ,j : 0 at r| -l 


u d* u 


du du 

*n ’ 0 a>,d " «>* ' Re (ft)* <*n J 


at -0 


(10a,b) 


( Ha , b ) 


where 5 - 6*/ b’ Q , H = u/u^, p - p/ x « x/6* 0 and Re * f ao u 0 o8o /, *ao 

Moreover, if wr tranefc rm x into 4 by the relationship 

x 


<U) 


o 
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Ihe last boundary condition becomes 
. dn Z d*u 

u 71 * \l) 1 77* at n o Hi) 

Equations (9), (IQa.b), (l la) and (13) are explicitly independent of 
p and p . Therefore, we can solve for u as a function of £, h fulfilling 
Equation (9) and boundary conditions (tGa, b)* (l l a) and (13) and also the 
initial data at £ - x - 9. The effect of and p enters the problem through 
the initial data* through the values of 6 Q and 0, and through the inverse 
transformation from the t> n plane to the physical plane of x, r. 

If the initial profile for u, which also fulfills the boundary conditions 
Equation^ (10a,b) and (l la), can be approximated by a third order poly¬ 
nomial in q as 

a* f j. n) - » - (* -*oH> - in*- 

1 1 4) 

or I - .1(0, r|) ‘ (1 - a 0 )(l-rj)*(l • ir\) 

then, the velocity profile in the wake can be expressed as 

I : [l-n(4>] ^-T ) ) (15) 

which obviously fulfills the conditions of Equations (10a,b) and (l la). 

Consider now the approximation of the initial velocity data of Figures 
1 and Z obtained with Equation (M). To make a comparison it is necessary 
to evaluate u(o, o) and 6 C * This was done by application of the method erf 
least squares. The resulting curves of l - u(o, r\) u m versus s are shows* 
in Figures 1 and 2. It will be noted that rough quantitative agreement is 
obtained, if a more exact representation of the initial profile is desired* 
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additional coefficients must be introduced into Equation (14); these must 
be treated as functions of £ with the requisite additional equations deter¬ 
mined either from additional boundary conditions at r * 0 or from the 
momentum equation multiplied b> u or r and integrated with respect to r. 
For the purposes of the present analysis the degree of approximation 
achieved by Equation (14) is considered satisfactory, Equations (9) and 
(i *) yield, respectively: 

;T •• 5* (o. IS(I-a).0.0857 (t-a)"*| 

° 8 i»6> 

~ o. i5(i-a o >-0. osafO-a^H* 

°o 

and a— - (U)~r- (17) 


Elimination of 6 yields 


~ * U(|.*) 4 fo. 15 - 0.085rfl-a3 
5 ’ d£ t™ -J 


By integration, £ ie expressed m terms of A or l-a(£). 


«o II 

£ r—r <0. 5556 (*T” - “ 0.^381 log 

6 o 1 . A *V> 


A 0 (1.75 -A Q ) 

A(|,75^o) 


where a 0 


Equations (i), (14) and (18) give the velocity and static enthalpy 
distributions in the £ - q plane, valid for either equilibrium or frozen 


flow and for either laminar or turbulent flow. These cases require 


separate treatment in the determination of the transformation from the £, q 
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plane to the physical plane x, r, of the transport properties and the electron 
density, ami of the initial values. 

Consider now the four cases cor responding to laminar and turbulent 
flow and to complete equilibrium and to completely frozen concentration. 

(A) F yiilibrium Laminar Flow 

Cor responding to each point i, t\ the static enthalpy h is obtained 
from Equation (Z) with the velocity u given by Equations (IB) and (14) * With 
p * the state of the fluid, the electron density n e$ and the viscosity 
are known for each value of h. 

At each station of x or 4* n e is maximum on the 4*axi«. The 
distribution of n e (^,0) with respect to x is obtained from the inverse 
transformation of Equation (!<;), 

4 

f' 

6 o R e J 

o 

which is carried out by numerical integration. 

(B) Fr ozen Lami nar F low 

In this limiting case, we assume that the volume recombination 
rates within the wake are low compared to diffusion across streamlines. 

The mixture in the wake in general is composed of molecule* O,, N*, NO; 
atoms O, N and electrons and the associated ionised molecules and atoms. 
The number of the ionized particles per cubic centimeter is large (cm the 
order of !0 J ^ for the conditions selected for numerical analysis), yet it ii 
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small compared tu the total number of partuses per cubic centimeter. 

Thus to a satisfactory degree of approximation the effect of ionized 
particles tan be neglected in the determination of the states of the mixture. 

If Y represents the mass fraction of species i in the wake, the 
continuity equation for this specie according to the integral method becomes 
the law of conservation of mass for ra* h species m the wake region, that is*,* 


, (ix r d u ^ ^ i to ~ ^ i 


) dr 0 


Since Y * 0 tor atoms‘and Y # 0 for mnlecules this- equation after 
\ OC l iX.i- 


transformation frurn r to* q becomes: 

# ,f 

f/t rj * \ dr, - l utj. t. 


► ub»K ript ; for atoms (*I0a) 


and 


(1 r \ q 

c 

l 

f£ f )*j r|;i (i ‘ Y k /Y ktJ> ) dq -const, subscript k for molecules 

UOb) 


o 


ihe bound.* tv Lonthtiuns ire 


cY 


at 


and at 


Y, 

n * o, 

dr dn 


i “ 0. *k'V| 
1 


>1 - I , Y - 0. - v . i k ' , k cw 

J d»| 


l, _>J5. , 0 

an 


la) 

Ulb) 


If Yj i» represented by a third order polynomial of q, the boundary 
conditions of Equations (ila.b) result in the dependence of Yj on q identical 


*lt should be noted that this law prevails independent of the relation between 
diffusions! velocities and concentration gradients. 
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to the dependence of I - u/u^ on tj, i.e. , 

b J<« u % 

Y, hji) (1-nV (l* 2b) —rrr (i - — ) 

J - i-aU) Uqp 


Equation UCa)then become* 

1 

b,(u c 

- L — \ n u( i - 

1 -aiU J 


— ) dr| = COfltft. 


From Equation { 9 ) it is clear that the above leads to 


b jts> _ 1 - aft) 

bjfO) 1-atOi 


u au~ a(4» b I 
Yj(O.O) " a ro - u(0,0) 


<22a* 


Similarly, if Y^/Y^^ is represented by * third order polynomial 
in Tj, the dependence uf l - Y^/Y^^ on rj will be identical to that of 

1 * «'»W ■ 


Yk C*> * Y j, {4 1 b) _ U qd- u (jL b) 
Ykai- Y k |0,0) u^-ufO.O) 


um 


It is of interest to reconsider at this point the assumptions with 
respect to the Prandtl and Lewis numbers. In the species conservation 


the assumption of a polynomial profile for Y* of the same degree a»|| 

velocity profile leads to Equation (22). However, if an additional coefficient 

^ f " 

' - 

resulting in a fifth degree polynomial lor Yj is assumed, and accor^B^f 

... * V-M* 

an additional boundary condition at r * 0 is imposed, there will result 
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differential equation analogous to Equation (17) for the determination of 
thin extra coefficient as a function of £• In this case the diffusion coefficient 
D t2 evaluated at the axis will arit*e* In an analogous fashion a stagnation 
enthalpy profile could be assumed as a substitute for the Crocco integral 
Equation (<£). The additional coefficient in this profile would be determined 
by a i or responding boundary condition at r ~ 0 arising from the energy 
equation. In this ase the Prandtl and Lewis numbers along the axis of the 
wake would ari^e, The initial values of the additional coefficients would be 
selected so that the initial profiles would be more closely approximated 
than is possible with the profiles used here* 

Knowing Y (£, n) and n L we can compute the density uf tne 

mixture y. The temperature of the mixture is determined such that the 
static enthalpy h of the mixture agrees with that from Equation U). 

Because of the close similarity between the transport properties and the 
atomic weights ot oxygen and nitrogen* the transport properties of the 
m;>ture in the wake will be calculated by the formula for a binary mixture 


of atoms and molecules according to the suggestion of Penner (Reference 9). 
The transport properties of binary mixture are obtained from the following 


formulas (cf. Reference 10, for example); 


>*1 l0 ' 5 ~~ 77 ^ r ;~ •*“ '- 2 

<r* c, rr.r 

i » 5 



GENERAL APPLIED SCIENCE LABORATORIES, INC. 





GENERAL APPLIED SCIENCE LABORATORIES, INC. 


M a t 


polynomi.il of tj, th« solution should bo similar to that for Yj, 1 ,#., 


r i )i' (0, 0) ag p • 'i(U 

*.-»<»■«- .— - —*— * * — "•* \fc 1 

\A0, 0>; (C, n) ~ u(C * u ; 

The inverse transformation from n to x f r plane can be 
accomplished by numerical integration since we know u (£, q) and 

^rnix. 


In Figure 3 there is shown the mass fraction of atoms given by 
Equation (ZZ) with 6*. determined from the velocity profile (t*f. figure ill), 
and with V (0,0) determined by the method of least squares. In Figure 4 
the ratio of the number of electrons per unit mass in the wake to the 
corresponding number at the axis a» given by Equation (Z7) again w ith h* 0 
determined from the velocity profile and n e (0, 0)f?{0, 0) determined by the 
method uf lea.ct squares* For both of these profile# the same comment* 
as were made previously u»r the velocity profiles are applicable* 
fC) Turhulent Wake ( Fc^uii ibrium or F runen Flow) 


In the i, n plane, the state in the wake is independent of p, 
therefore, the result# for equilibrium or frozen laminar wake flow* can 
be employed for the corresponding turbulent flow provided that the symbol 
p is interpreted as the absolute eddy viscosity p f instead of the absolute 
viscosity p in the transformation from n plane to the physical x-r plane. 
Experimental data for absolute eddy \iscotity for compressible wake flows 
du not appear to be available, Mager {Reference 11) has shown, however, 
that for boundary layer flows the compressible values of p* can be 
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correlated ti> the in'ompressibl* eddv vi^co ity **. Mager's correlation 
can be mothfi'-d t’sec* Appendix) for the transformation of compressible 

I 

axisymmetric wake flow; there result 

r 1 . j 

^ , i mm 1 

P t - ~—(, P ci r i (48) i 

9 1 J 1 

o . ■ - - ' i 

According to Reference l, pages 4bd-S01, the incompressible eddy viscosity 
* 

«* for wakes and jets can he expressed as 

£w u<0 '°2) w\ 

where (6^)^ : 6^ is the radius of the cor responding incompressible wake 
at the initial station. 

The value of ranges from 0. i)£ tv 0.005 for mixing of two* 
dimensional or circular jets and for two-dimensional wakes* Since a 
different value of k amounts only to a different seal* of k in the physical 
plane, the value of 0.0! considered typical will he employed in the present 
calculation. With p replaced by p t the procedure developed for the laminar 


wake can be applied tv.* the turbulent wake. 
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IV ELECTRON CONCENTR ATI ON I X A LA MINAR, SHOC K -IN DUCPD 

WA KE 

The traii-'l«»rmatun.s to the physical plane h«‘tve presently been completed 
tor the cases corresponding to cases A and B oi the previous section* The 
results in terms oi number of electron* per unit volume along the axis 
n e (o # x) are shown it Figure 5. It should bo noted that the x scale therein 
corresponds to a hasp diameter ot one too*; for a base diameter equal to 
2.5 feet these dimension* should be multiplied by (2.5)* » o. 25. 

For the* case oi equilibrium flow the results of Figure Syndicate that 
the number of electron© per unit volume relatively low at the initial 
station and decrease® by a factor of two in roughly 24,000 feet lot a body 
of base diameter 2. 5 feet. For frozen ionization the number of electrons 
per unit volume is relatively high at the initial station U.8 x 10 tl ) and 
decreases by a factor of two in 9, 200 feet, again for a base diameter of 
2*5 feet* 
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\ FFECTR<ON a >NC FN I P A UOHltl A rUKHULENT, SHOCK - _ 

INDUCED WAKK 

Rased upon tut* trundormatiuafi tx* Section III C, the results m turmi 
ui number at elections per unit volume along the axih n e (m x) are shown 
in Fit; .re t for b«*th equilibrium flow «nd frozen flow. 

As l umpar n to the < or respond; ng laminar cases m Figure S f the 
length for the ><xti\e d« c*iy in electron density is redrn id by an order 


of Af>0 fur the equilibrium flow and by an order of dbO for the Iroattm flow 
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I he anaiysi. for the case of boundary layer induced wake can be 
carried out along lines similar to that described above for the shock* 
induced wake. Wo consider as a first analyses the approximation of 
Prandtl number of unity. Thet/since the wake is treated as a constant 
pressure region, the Crocto integral given by Equation (2) can be extended * 
so that 


h * u^/2 » A ♦ B u 





whore A and B are constant* de termine d by : 
clearly that h * h m when u * u<# whti* the-' : ^/' s ■ *'" /‘ : j 
a specification of the total energy lu*u by the stream to the body generating 
tn<; woke. Consider* therefore* that A and B are known. 

the solution for the velocity distribution on the qi plane .given 

Equatiuna (U) and il$) is still valid for the case. rtowevetv the etrthfipy 
distribution will now be obtained from Equation O0) rather than from 
Equation U) me therefore will change* The transformation to the x, t 
plane obtained from the inverse of Equation (12) must also be carried out 
numerically and wail be altered by the change in the velocity*<mthhJ|py 


relation. 
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vu numer i cal exam ple of boundary layer induced wake 

In order to provide an indication of the wake associated with the 
boundary layer, the flow about a slender, spherically-capped cone is 
considered. The half-cone angle is taken to be iS° so that the ionisation 
and dissociation behind the conical shock is negligible. Moreover, a 
base radius of £* S feet is assumed in order to permit comparison of this 
boundary layer induced wake with those computed above. Finally, the 
radius of the spherical cap is selected ho that at the flight condition# cor* 


responding to £00*000 feet the mass entering the boundary layer at base 
of the body has passed through the . conical shock, Theboundary 
conservatively assumed to be laminar.# 

In calculating the radius of the spherical cap the Levy-Lees tra»s~ 
for mat, on (cf* Reference 10) is applied along with the assumption of uniform 
external flow on the cone (0 W OK The mass entering the boundary lay## 
at a station x measured along the surface of the body originates m a 
stream tube of radius according to the following equation 


P m V <x>* s ^ * N 2# lfh e > V- ; 

where the subscript m denotes flight conditions, y i* the stream. 

s’is the transformed x coordinate defined by the equation 

‘ • ‘ 

V • ■ ' . , , ‘. ^ . .?$&• *i£h* : 
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and finally where f( i| ) is the value of the transformed stream-function 

evaluated at the "edge 1 * of the boundary layer. 1 * For u/n^ fi.994, f<h*) *? £♦ ** 

It is sufficiently accurate for present purposes to neglect m the rumpu* 

^ I il 

t&tion of s the influence of the spherical tap which is exi opted *o be ot 
relatively small radius. Thus using a Newtonian pressure nn the cone and 
the oblique shock properties given by Reference 4, we find from 
Equation; (31). There was available a shadowgraph obtained at a Mach 
number of 7.9 on a spherically-capped I :>° half-angle cone. Employing 
Mach number independence principle, we can find the value of the no 
rauius such that the sireamtube with the Computed value of t m * {r^j^ 
passes through the junction of the conical and curved portions of the bow 
shock. This nose radius is found to be 4. i inches* 


With the nose radius determined, the shadowgraph can be used *n 
conjunction with the gas and flow tables of References 4 and 5 to determine 
the state and composition of the gas behind the shock and thus as a function 
of the radial position of the point in question on the shock* 1 he mass 
fraction of the constituents considered are shown in Figure 7. The**** data 

are required for the determination of the initial condition in the wake 

analysis* 

Consider next the thermodynamic process experienced by the gas {$* 
flowing from the region behind the shock, through the boundary layer a&d 
into the wake behind the body v The process leading to the greatest electron 









population in the wake ns as in the cast- of the shock •induced boundary 

that of completely froaen flow; thus for the subsequent analyst* thiedahrer* ' 

vativ*; case will be considered. 

We discuss first the miti.il profiles of velocity, density,,, 
and composition in the wake where as above we take p a; p^v, ; 
and stagnation.enthalpy distributions at the trailing edgey>i 
assumed to be given by flow on the cone .<»! * 3 ); in -|fr«i*a 
transformations with the assumption* pp*»•«,**;» n*Kt:hai^.;-F.^ 
the Blasius values tabulated 


■employed* - With the ass«tmp*i#n; : df-‘ : id 




edge of the cone to the pressure p , the velocity u(o, of at ^e dd»tdrr".ol-;^fet:- : 
initial cross section of the- wake is obtained. At the edge of the wake 


r =6 u or it • ^ u « where ^ is the,associated, 
radius as defined'in Equation fSel. Witb the 'ftji■>|rinijf" 
the velocity profile equivalent to Equation (IS), i.e., 


ufo, q > - u(o, o) 
u(o, 1) « u{o, o) 


* (I - q) 4 (I + <iq) 


where q » the initial wake radius 6£ is determined by the condition-«|g 


conservation of mass: 


p «o -\ ufo,q)qdq • \ purdr * 


mmm 


It will he noted that the initial radius Is l. 04 feet ** 
in the caee of.the shock l»duced^wike fe»';i:fhp^M,i*:lf ; '' 


: ,r x :>/• 

yj'5 tu ■’‘•Si’J'v 1/ r-'rsf? 


I 




-c- 1 








- 'rl’lTr'i 
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The farther determination of the g <-& composition is more complex. 
Consideration of Uw* *nass fractions of constituents shown jn Figure ? 
indicates that no Kim pie poh.i omul expressions c*.n he employed to 
approximate the relations between mass flow and masa fraction, for the 
purpose of estimating the initial state of the wake, straight-line segments 
were employed along with the requirements of species conservation 
corresponding to frozen flow; this requirement can be expressed as 


U u 
' f uo ou 


Y ; u r dr) - f V p ¥; u r dr) 

shock J -x*o in wake' 

O...... • .. , 


The left hand side of Equation (52) can be evaluated for each i from 
Figure 7 ? while for the riant h^nd .% x de it is convenient to use the traus- 
farmed variable q in place of r, thus from Equation (8a) 


{ » r V.- u r cr) 


x*o m w ake 




x * a its wake 


Provided straight line .-ogments for Yj {q ) are assumed the distributions 
of the mass fractions of the main constituents shown in Figure 8 are obtained. 

We are now able to find the density and temperature distribatione 
car responding to the pressure Pqq and to the static enthalpy distribution 
which ia given by the Croeco relationship since -r.;* T 


p - uu - 


R T(t) 


Vhi 
L. w 


Y i(U 






GENERAL 


* hit* ‘I {.-) it- given pafSriSeiri^ny’1>y theecquati on 


h * h<? * S X h ‘ ( 15 


The distribution oi pfp^ v** 1 ! is shown In Figure 9. 


In * rdf r tu e.-timaie the electron distribution, it wilt be app r ox i nitt(R Ly 


a third order polynomial of q, i, e., 


n JO, »)}* *dO,m ♦ in) 


due to the boundary conditions at if*' 1 and q ~ G. 


t he «m;l»tude %(6,.OMN^-nntwM;by tte condition of conservation 


total electron populationfrom the shock tothewake. 






\ * j.• *,v 'v ;. fti: ,-c*>’V if- ijyjfi&zi 


n e r i» <i? 00 * < V* 4 n e W * %M p 1 * dT 2 ■ 


This condition yields a higher estimate because of the loss of electrons La 


the stream tubes near the surface of the cone. 


of that for a shock induced wake. The initial radius i« approximately 1/9,4 


times that in the shock induced wake. The decay of electron density can be 


estimated front that in a shock induced wake wi th the length scale along the 


ilsi 


d*c r« a » e (torn I«II :x 


! 


: V"\ 


eM 


\ 


■ 
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VIII CONCLUDING REMARKS 


An analysis of th*. fluid mechatmai aspects of the wake* behind bodies 
in hypersonic flow has been carried ovt* The fluw has been idealised so 
that the pressure is constant and the boundary layer approximations are 
valid. Both laminar and turbulent flows aie considered. Furthermore, 
two limiting cases are considered corresponding to the shock-induced wakt 
and to the boundary layer induced wake. The main interest in this study :> 
the electron concentration in the wake. Preliminary analysts indicated that 
the greatest number of «Uctrw$ would arise from tb# ca#e# of complete 
equilibrium flow or completely frozen flown therefore, the chemical and 
thermodynamic behavior of the gas is idealised m two different ways 
corresponding to these two cases* 

By appropriate transformations the velocity distribution m a trans¬ 
formed plane can be obtained for all cases of laminar or turbulent flow, 
shock- or bouncary layer-in luced wake, frozen or equilibrium gas behavior. 
The inverse transformation to the physical plane differs m each case and 
must be performed numerically# 

Numerical results have been obtained for lamina? and turbulent flow 
with a shock “induced wake corresponding to the 3* IQ configuration. The 
flight conditions were taken to be £3,200 fpa and 200, Q00 feet altitude. 
Although the detailed calculations were carried out for a base radius of ono 
foot, the scaling to a base radio# of l, 5 feet is easily performed, 

Of equilibrium, laminar flow t# found to yield an ettctron conce the 


mM 




mtmAL Amm s cience iA» 0 M 70 tits r me. 


%z 


origin at the wakv of roughly d x 10* .e!*ctron*/cc# This concentration 
decays by a factor of two m a length of 44, 000 feet for a body with base 
diameter of 4. 5 feet* For fro?«ri flow electron concentration is 

j 

approximately 1,8 x iu *7cc at the origin uf the wake and decays to one*half 
in a length of b, 000 feet for a body of the same base radius* For ih# 
corresponding turbulent cases, the scale of length is reduced from the laminet 

values by roughly >50 fur equihbirum flow <md by 450 for freren flow* - 

. ■ ' ■ . . ' , v / I : ; " a .1 , ' 

As an indicate?) oi the wake associated with the boundary layer, the flow 


about a slender *»pheru ally > apped cone i-%-_cpr t '*ide.r^d'*-. ^3!&* tp&f. 

> • • • -'•••. '■ 1 . ' ' *•" ' f <■ i'-V^, yt y-rt'‘ 'C 

is15°* the base- radius t-e 4/5 i e«t;and.• thecae* 

chosen so that the fluid passing through the curved position of the bow shock 
just com pi ♦•tod is -w allowed b\ the boundary ’aver at the trailing edge of the 
cone. The maximum electron density.at the initial crossrsec.tfo^; of-. y-i 
wake is -roughly bi x 10 ptrc*c v . It i* one-fifteenth of 

value in the shock-induced wake corresponding to the iyitt-'. : y : 
decavs to one-hali m a distance of 100 feet. 
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a pprr:r; y * 

i KANSF ORMAT ION OF i HE T URBULENT 
COM P R ES S1HL E WAKE FI OW 

1 he g r ^ermng ^quatians for axi symmetric turbulent wake ficiw sire: 
(pur) x t (^vr 4 p f v*rj y ~ 0 (A-1) 

rpu u x + {pv t * p , v t r)u r Mpr u f v f ) r 

purH x * (pvr t r)H r - r <A-3> 

We hr*v *» already employed two assumptions that the pressure is nearly 
con- tant m the wake and the effects of viscosity 
ner^h‘1 t.;d 

In these equations the primes rienntc the fluctuating terms and H # th# 
total enthalpy, includes the mcar kinetic energy* Equation (A-3} t» fulfilled 
by th»' solution H« constant and v’H* * 0. 

In addition to the requirement m the laminar flow transformation Att 
the stream function remains unchanged, we impose on the transformation 
the following condition*: 

a. The reference condition* at the free stream value should be need. 

i>, The moment of the turbulent shear about the axis is constant over 
an infinitesimal mass of the fluid. 

Condition (a) is more appropriate for wake flow than Mager*s which 
was introduced for boundary layer flow. 
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Condition (b) ib aldifferent from Mager** because ol the appearance 

of the moment of htrens rp uH * in the equation of motion* 

* 

With ( ) indicating the quantities in the transformed > vstem the 


transformations are 


(A - 4) 


r 

(r*)* * j (JL ) dr' 2 * z* 

t, P«, 


(A-5) 


(A-b) 


4» - 


f A - 7} 


(r,' u*\ ’) p (rdrdx) * it a u'v* ’*) ,, r*dr*dx* 
* c© OO 

The la?i equation lead# to 


(A-f) 


u*v f r r**> u^v* * 


{ A - 9) 


The stream functions art- related to the velocity component:-: a$s follows: 
pur 

and u*r* * >i»p* v*r* - -4»^* 

Equation (A-7) yields 

u* - u and r{<iv + p’v 1 ) * i'^ |r*v* - r*r^u*J 
The continuity and momentum equations for the compressible wake flow 
(Equations l and 2) then become the corresponding equations for incompres¬ 


sible flow: 


(u*r*J x *+ (v*r*) r * * 0 


(A-10) 


u * u ** + v * u *# * 'JTi (r* u'v* *) 


(A*U) 
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If the eddy vi He osi tv i :> introduced, F equation (A *9) becomes 


d.i 

r in *' - r* p f *- 

r or uo ax* 


or 


r* 2 Kciu Pun* * 

~~r" 

f p r j 

o 


P d r* 


Along the axis, we have p Jx, o) * p <’ 
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Radial Distribution of Mans Fraction of Major Constituents 

Behind Dew Shock 
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